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Abstract. Precise radial-velocity observations at Haute-Provence Observatory (OHP, France) with the ELODIE 
echelle spectrograph have been undertaken since 1994. In addition to several discoveries described elsewhere, 
including and following that of 51 Peg b, they reveal new sub-stellar companions with essentially moderate to 
long periods. We report here about such companions orbiting five solar-type stars (HD 8574, HD 23596, HD 33636, 
HD 50554, HD 106252) and one sub-giant star (HD 190228). The companion of HD 8574 has an intermediate period 
of 227.55 days and a semi-major axis of 0.77 AU. All other companions have long periods, exceeding 3 years, and 
consequently their semi-major axes are around or above 2 AU. The detected companions have minimum masses 
7712 sini ranging from slightly more than 2 Mjup to 10.6 Mjup. These additional objects reinforce the conclusion that 
most planetary companions have masses lower than 5 Mjup but with a tail of the mass distribution going up above 
15 Mjup. The orbits are all eccentric and 4 out of 6 have an eccentricity of the order of 0.5. Four stars exhibit 
solar metallicity, one is metal-rich and one metal-poor. With 6 new extra-solar planet candidates discovered, 
increasing their total known to-date number to 115, the ELODIE Planet Search Survey yield is currently 18. 
We emphasize that 3 out of the 6 companions could in principle be resolved by diffraction-limited imaging on 
8m-class telescopes depending on the achievable contrast, and therefore be primary targets for first attempts of 
extra-solar planet direct imaging. 

Key words. Techniques: radial velocities - Binaries: spectroscopic - Stars: brown dwarfs - Stars: extra-solar 
planets - Planetary systems - Stars: individual: HD 8574; HD 13507; HD 23596; HD 33636; HD 50554; HD 106252; 
HD 190228 



1. Introduction lar companions and allowed them to discover the first 

' extra-solar planet (Sl Pegb) orbiting a solar- type star 

. - . The ELODIE Planet Search Survey, an extensive radial- jMavor fc QuelodliS . This discovery was followed by 

velocity northern survey of dwarf stars, has been under- ^ ^^^^^^er of subsequent detections of extra-solar planets 

way for several years at the Haute-Provence Observatory ^-^^ ELODIE. Some of these results have previously been 

(OHP, CNRS, France) using the ELODIE high-p r ecision published 

fiber-fed echelle spectrograph ( |Baranne et alj E^M) This survey is part of a large effort aiming at 

mounted on the Cassegrain focus of the 1.93-m tele- gxtra-solar planet search through radial-velocity mea- 

scope. This survey was initiated in 1994 by M. Mayor g^rements since several other planet surveys with a 

and D. Queloz in order to detect very low-mass stel- sensitivity Ki>10ms-^ are underway elsewhere. A 

list of the most productive ones includes the Geneva 



Send offprint requests to: Christian Perrier, e-mail: Southern Planet Search Programme, carried out with 
christian.perrier@obs.ujf-grenoble.fr CORALIE spectro graph mounted on the 1.2-m 

Based on observations made at the Haute-Provence ' ' ' ' ' 



Euler Swiss telescope jOueloz et al.l l2000bt lUdrv et all 



Observatory (operated by French CNRS), the 1.2-m Euler U^I I , x^n^ x n-n^ .i ^ '. 

. , , , T-,r,r^ T c-11 r^u 4- ^r^u-i \ j I2000al l at ESO-La Silla observatory, the California & 

SWISS telescope at ESO-La Silla Observatory (Chile) and the ' ^ 

1.52-m ESO telescope also at La SiUa Observatory. Carnegie Planet Search with the HAMILTON spectro- 

The ELODIE measurements discussed in this paper are S^aph at the Lick Observatory (|Marcy fc Butler, ^993) 

only available in electronic form at the CDS via anony- and the HIRES spectrograph mo unted on the 10-m 

mous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via Keck-1 telescope (|Vogt et al .l ll994j) at the W. M. Keck 

^http : //cdsweb . u-strasbg . f r/Abstract . litml_ Observatory, the G-Dwarf Planet Search also performed 
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Table 1. The 18 planet candidates with m2sini < 
11 Mjup detected with ELODIE at Haute-Provence 
Observatory 



Candidate 


References 


51 Peg b 


Mavor & Oueloz (1995) 


14 Herb 


Udrv et al. ('2000b~l. Naef et al. f 2003a) 


Gl 876 b* 


Delfosse et al. (1998") 


HD 209458 b* 


Mazeh_etaL (2QQU) 


HD 190228 b* 


Sivan et al. (2000) and this Daoer 


HD 8574 b* 


this paper 


HD 50554 b*"^ 


this paper 


HD 74156 b* 


Naef et al. f2003a) 


HD 74156 c* 


Naef et al. (2003a) 


nu OUDUD u 


i > aei cL ai. izuuiai 


HD 106252 b*t 


this paper 


HD 178911 Bb°* 


Zucker et al. (2002) 


HD 20367 b* 


Udrv et al. (2002) 


HD 23596 b* 


this paper 


HD 33636 b*t 


this paper 


HD 37124 0*" 


Udrv et al. (2002) 


HD 150706 b* 


Udrv et al. (2002) 


GJ 777 Ab*^ 


Udrv et al. (2002): Naef et al. (2003b) 



See also 
See also 



Marcv et all Jl99. 



Charbonneau et al.l (|200(t : iHenrv et alJ ||2Q0 



* Announced in ESQ Pre ss Release nb 07/01 (April 4, 2001) 
^ Confirmed bv lFischer ct al. ( 2002) 

° Variability detected by the K eck G-Dwarf Plan et Search 

* Independently announced by IVoet et all (|2002h 

* Announced during the Scientific Frontiers in Research on 
extra-solar Planets conference, Washington, June 2002 

See also Bu tler et al.. (.200a) 
^ HD 190360 b 



with the HIR ES spectrograph at Keck Observatory 
l)Lathamll200(]() . the Anglo- Australian Planet Search with 
the UCLES echelle spectrograp h mounted on the 3.92- 
m Anglo- Australian T elescope llTinn^ y e^ ^iJ I2OOII) . the 
AFOE Planet Search teorzennik et all 119981) using the 
AFOE spectrograph mounted on the 1.5-m telescope at 
the Whi pple Observatory and th e McDonald Planetary 
Search ijcochran fc HatzesI Il994j) using the McDonald 
Observatory 2.7-m telescope coude spectrograph. 

115 planet candidates, including 11 planetary sys- 
tems, with minimum mass es lower t han 18 Mjud and down 
to su b-Saturnian planets ijMarcv et al.. .2000; Pepe et alJ 
120021) . have been detected so far, with periods from slightly 
less than 3 days to about 15 years. This should permit us 
to obtain more robust conclusions about the mass function 
estimation, and notably its fast rise towards lower masses, 
the relative metal content of the parent stars and con- 
straints on t he so-called brown dwarf desert at separations 
of a few AU ijHalbwachs et aLlbflOfll) a nd on the ernergin g 
period desert between 10 and 70 days 1 Udrv et alJ l2002V 
crucial ingredients to the understanding of planet forma- 
tion mechanisms. 

This paper is the first of a series devoted to the dis- 
covery of extra-solar planet candidates in the framework 




6000 



5500 



5000 
, (K) 



4500 



Fig. 1. The Hertzsprung-Russel diagram for the ELODIE 
main programme sample stars. The six stars with new 
planet candidates are shown as open triangles. The four 
curves are the solar metallicity evolutionarv tracks for 0.8, 
0.9, 1 .0 and 1.25 Mq from Geneva models I Schaller et alJ 



of the ELODIE Planet Search Survey. It reports on the 
detection of six new candidates, with minimum planetary 
masses (m2 sini) in the range 2 to 11 Jovian masses. Five 
of these are orbiting F8 or GO dwarf stars and one is found 
around a G5IV star. They are all but one long-period ob- 
jects on little to moderately elongated orbits: three periods 
are found to be greater than 1500 days, the shortest one is 
227.55 days and the orbital eccentricities are of the order 
of 0.5 for four objects and close to 0.3 for the two other 
ones. These discoveries, together with 3 others announced 
during the Scientihc Frontiers in Research on extra-solar 
Planets conference in Washington (June 2002), bring the 
number of planets discovered with ELODIE up to 18. 

The next section describes the sample of the ELODIE 
Planet Search Survey and summarizes the planet candi- 
dates it has allowed to detect up to now. 

The basic properties of the stars harboring the six 
new companions are given in Sect. 13 Sect. ^presents the 
radial-velocity data and the orbital solutions for the com- 
panions. The case of HD 13507, suspected for some time 
to harbor a planet but that actually has a low mass star 
companion, is exposed in Appendix A. 

2. The ELODIE Planet Search Survey 

The original sample of iMavor fc Ouelod l)l995|) contained 
142 stars, out of which 51 Peg. The main selection crite- 
rion of these G and K dwarfs was their radial-velocity 
non-variability according to the anterior CORAVEL sur- 



C. Perrier et al.: The ELODIE survey for northern extra-solar planets I 



3 




150 



100 



50 



— 1 1 1 1 — : 



_I I I I l_ 



0.5 



1 1.5 

B-V 



Fig. 2. The histograms of my and B — V for the ELODIE 
main programme sample stars 



vey i Duaiiennov et al]ll99lt IPuanennov & Mavo3ll99ll) . 

The sample was largely modified in 1996 and later, so 
that the to-date survey sample size amounts up to 372 
stars. Among the present sample, only 71 stars remain 
from the original list (active and/or too faint stars of this 
list were excluded) and 259 additional stars were included 
with the following criteria: my < 7.65, (5 > °, v sin i (from 
CORAVEL) <4 kms^i, CORAVEL non-variability, spec- 
tral type from F8 to MO. Both lists are now forming the 
main programme with 330 objects. 

In addition to these main programme stars, for follow- 
up or phase-coverage reasons based on the large access one 
can have to a 2 m-class telescope observing time, the sam- 
ple was more recently complemented with 42 stars com- 
mon with other programmes, among which 32 stars, with 
V magnitudes up to 9, were id entified as vari able by the 
Keck G~Dwarf Planet Search (|Lathamll2000|) . One star, 
HD 209458, was in our main programme and now also be- 
longs to the collaborative effort with the G-Dwarf Planet 
Search. 

The Geneva Southern Planet Search programme and 
the ELODIE Planet Search Survey have distinct stars sam- 
ples, based on the declination, although the former may 
occasionally support the monitoring of a usually faint, low 
declination ELODIE programme star. 

The planetary companion detection is based on 
the stellar radial-velocity reflex motion produced by 
gravitational interaction between companion and star. 
The ELODIE spectrograph, operating between 3850 and 
6800 A, has a spectral resolution of about 42 000. The 
planet survey makes use of its dual-fiber mode where a 
star and a reference thorium-argon calibration lamp are 
simultaneously observed to monitor and correct the in- 
strumental drift and to achieve the required high preci- 
sion. The signal-to-noise ratio per pixel for my = 7.65 
stars is about 100 for a 15-minute observation under good 
atmospheric conditions at OHP. In order to achieve the 
most homogeneous observations, some exposures are du- 
plicated in order to reach the target signal-to-noise. None, 
except for some special cases, are longer than 15 min- 
utes in order that no differential instrumental drift be- 
tween the two beams dominate the resultant accuracy. 
The total exposure time per source is dictated by the 
goal that the photon-noise be lower than the instrumental 



limit. The radial- velocity measurement is produced on- 
line by the ELODIE automatic reduction software which 
cross-correlates the observed stellar spec trum with a bi- 
nary mask template ( Baranne et al.ll996l) . Some improve- 
ments have been implemented in 2001 in those regions 
of the mask partic ularly rich in tell uric absorption fea- 
tures (see for details lPeoe et al.ll2002j) leading to a current 
instrument-limited precision better than 6ms^^. 

The improved precision is demonstrated by the dis- 
tribution of the weighted radial- velocity rms measured for 
111 constant stars in our sample (Fig.lSJ. We corrected the 
measured dispersions for the contribution of the photon 
noise. The instrumental error can be well estimated by the 
rms values of the low-dispersion tail of the histogram. The 
median of this distribution is 7.8 ms~^ and we find that 
20% of these constant stars exhibit a velocity dispersion 
lower than 6ms^^. We also show in Fig. 0] our data for 4 
constant stars with dispersions lower than 6.5 ms~^ (pho- 
ton noise included). Finally, the residuals to the orbital 
solutio n obtained with ELODIE for 70Vir ijNaef et alJ 
l2003a() are 6.1 ms~^ (photon noise included). The photon 
noise error of the 70Vir measurements is ephot = 2.8ms^^. 
The instrumental error inferred from this long series of 
values is: Cinstr = (6.1^ - 2.82)i/2 = 5.4ms-^ 

More recent improvements in the telluric absorption 
features treatment have allowed us to reduce systematic 
yearly velocity effects that were still visible for some stars 
in our sample. This is for example the case for HD 33636 
where the residuals to the fitted orbit could be reduced by 
a factor of almost two. 

A total of more than 9000 measurements, spread over 
observing periods that are programmed about once every 
month, have been performed since the beginning of the 
programme at OHP. All stars but a few ones have been 
observed more than 10 times and in average almost 20 
times but those presenting indications of variations may 
have been observed more than 50 times. 

Table ^ lists all the planet candidates discovered so 
far with ELODIE having r7i2sini< llMjup including the 
six presented in this paper and four others also recently 
announced. More massive companions, with minimum 
masses in the range 20-75 Mjup, among which HD 127506, 
HD 17445 7 and HD 185414 , have been discovered with 
ELODIE l|Naef et al.ll2000bl Naef et al. in prep). A pecu- 
liar case, HD 166435, has also emerged from the ELODIE 
Planet Search Survey observations. This star presents a 
chromospheric activity stable enough to mimic a short 
period Keplerian signature. Therefore, any radial-velocity 
variation potentially due to a short period companion is 
now first checked against such activity using line bisector 
va riations, tha t have been shown to be a robust diagnostic 
bv lQueloz eral. (.2001.) . 

The hst in Table [H includes HD 209458 b l)Mazeh et alJ 

l200nl) . the only extra-solar planet for which tran- 
sits we re observed. Photometr ic transits were ob- 
tained ijCharbonneau et al.l l2000() with timing predic- 
tions computed from ELODIE monitoring and commu- 
nicated to D. Charbonneau and collaborators. The first 
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Fig. 3. Histogram of ELODIE radial-velocity weighted 
rms for constant stars. The 1723 measurements of 111 
stars having P(x^) > 0.05 and iVmeas > 5 have been used 
for the plot. These selection criteria are sufficient for elim- 
inating highly active stars but lower activity objects can 
still be present in the sample. The solid curve represents 
the distribution of the weighted rms corrected (quadrat- 
ically) for the photon noise. The main error sources still 
present in these corrected dispersions are the instrumen- 
tal error and the stellar jitter. The median of this distri- 
bution is 7.8 ms~^ and 20% of the stars in this sample 
have weighted rms lower than 6ms^^. This latter value 
value is a good estimate of the ELODIE instrumental error. 
We show for comparison the distribution obtained without 
correcting for the photon noise (dashed line). 



spectro-photometric transit was observed with ELODIE 
( Queloz et al.ll200o'al) . The photometric and spectroscopic 
transits yielded notably the planet radius, the rotation di- 
rection with regard to that of its primary and the mean 
density (0.3 gcm""^), essential to assess the gaseous nature 
of the planets. Two planet candidates in Table were 
found around stars not belonging to the ELODIE main 
programme: HD 178 911 Bb, a planetary companion in a 
stellar triple system (jZucker et al.ll2002|) . and HD 80606 b, 
the extra-solar planet having the highest orbital eccentric- 
ity known so far ( Naef et aL 2flQla,) . They are both belong- 
ing to the G-Dwarf Planet Search programme carried out 
at the Keck telescope. Also, G1876b, which is presently 
in the closest extra-solar planetary system known, was 
discovered as a by-product of an M-dwarf multiplicity 
statistical study carried out with the OHP-ELODIE and 
ESO-FEROS spectrographs and with the Canada-France- 
Hawaii Telescope (CFHT) PUE'O ad aptive optics imagin g 
system by a Grenoble -Geneva team ("Pelfos se et al Jl99^, 
and independently bv lMarcv et al.l (,1998..) . Later, the pres- 
ence of two planets on resonant orbits with periods in the 
ratio 2:1 that were masquera ded by a un i que, m ore eccen- 
tric one, was recognized bv .Marcv et al ] il2nnih . 
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Fig. 4. Radial- velocity measurements in kms""'^ versus 
time for four standard stars of the ELODIE programme 
sample: HD 18803, HD 97711, HD 109358 and HD 117043. 
The mean photon noise (ep„) is indicated at the lower 
left of each frame and the measurements photon noise- 
weighted rms error at the lower right. The time span of 
these measurements is similar to that for the six stars with 
planet candidates. 



Fig. n] shows the Hertzsprung-Russell diagram of the 
317 stars in the ELODIE Planet Search Survey main pro- 
gramme sample for which an HIPPARCOS entry exists. 
The six stars harboring planet candidates discovered with 
ELODIE, presented in this paper, are shown with distinc- 
tive symbols. Tcff is derive d from the HIP PARCOS B-V 



by using the calibration of iFlowe^ (^2^^, except for the 
six planet-candidate hosts for which the effective tem- 
peratures listed in Table [3 have been used, and L from 



^ Note tha t the coefficients of the calibration printed in 
iFloweJ (Il99d) are incorrect; we rather used the coefficients di- 
rectly obtained from the author. 
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HD 8574 HD 23596 HD 33636 



HIP 




6643 


17747 


24205 


Sp. Type 




F8 


F8 


GO 






7.12 


7.25 


7.00 


13 — V 




U.0< / ± U.Uii 


u.Do4 zh u.uuy 


U.ooo ± U.UiD 


TT 


(mas) 


22.65 ± 0.82 


19.24 ± 0.85 


34.85 ± 1.33 


Distance 


(pc) 


44.2 ± li 


52.0 ± y 


28.7 ± 1.1 


/Xq cos(o) 


(mas yr ) 


252.59 ± 0.76 


53.56 ± 0.68 


180.83 ± 1.07 




(mas yr ) 


— 158.59 ± 0.55 


r\r' 1 r\ K 

21.06 ± 0.55 


— 137.32 ± 0.64 






3.89 


3.67 


4.71 


n.C 




n no ,1 
— U.Uo4 


— U.U29 


-0.046 


L 


(L©) 


2.25 


2.76 


1.07 




(K) 


6080 ± 50 


6125 ± 50 


5990 ± 50 




(cgs) 


4.41 ± 0.15 


4.29 ± 0.15 


4.68 ± 0.15 




(kms-i) 


1.25 ± 0.10 


1.32 ± 0.10 


1.22 ± 0.10 


[Fe/H] 




0.05 ± 0.07 


0.32 ± 0.05 


-0.05 ± 0.07 


Wx,u 


(mA) 


52.1 


76.7 


49.5 


log n(Li) 




2.56 ± 0.09 


2.81 ± o;j7 


2.46 ± 0;;^ 


V sin i 


(kms-^) 


4.04 ± 0.61 


3.59 ± 0.59 


2.79 ± 0.65 


M. 


(Mq) 


1.17 


1.30 


1.12 






\\ \ ) OU0O4 






HIP 




33212 


59610 


98714 


Sp. Type 




F8 


GO 


G5IV 


mv 




6.84 


7.41 


7.30 


R T/ 

_D — K 




U.Ooz ± U.UUo 


U.uoo =h U.UU/ 


n 70Q -1- n f\c\(\ 
U. ( yo =h U.UUd 


TT 


(mas) 


32.23 ± 1.01 


26.71 ± 0.91 


16.10 ± 0.81 


Distance 


(pc) 


31.0 ± 


O '~7 4 1 1.3 

37.4 ± 1.2 


62.1 ± Ig 


Ha cos((S) 


(mas yr ) 


-37.29 ± 0.93 


23.77 ± 0.91 


104.91 ± 0.37 




(mas yr ) 


—96.36 ± 0.50 


—279.41 ± 0.53 


—69.85 ± 0.57 


Mv 




4.38 


4.54 


3.33 


ij.G. 




—0.038 


—0.061 


—0.176 


L 


/T \ 

(Lq) 


1.45 


1.27 


4.31 


erf 


(K) 


fiO'iO + ^0 


"iSQO + 'iO 


"i^fiO + 40 


logff 


(cgs) 


4.59 ± 0.15 


4.40 ± 0.15 


4.02 ± 0.10 




(kms-i) 


1.19 ± 0.10 


1.06 ± 0.1 


1.12 ± 0.08 


[Fe/H] 




0.02 ± 0.06 


-0.01 ± 0.07 


-0.24 ± 0.06 




(mA) 


46.6 


9.7 


8.8 


log n(Li) 




2.48 ± 0.10 


1.62 ± gj? 


1.05 ± III 


V sin i 


(kms^^) 


3.32 ± 0.59 


1.74 ± 0.25 


<1 




(Mq) 


1.11 


1.02 


0.83 



HIPPARCOS apparent magnitude and parallax and the 
bolometric correction given by Flower's calibration. 

In Fig. 13 the histograms of my and B — V for the same 
stars as in Fig. ^ show a clear cut-off in magnitude dic- 
tated by the ELODIE performances and a range in spectral 
type with a pronounced bias toward G0-G5 dwarfs and a 
strict cut-off at the F8 spectral type. 

3. Stellar properties 

The six stars found to harbour planetary companions are 
HD8574, HD 23596, HD 33636, HD 50554, HD 106252 and 
HD 190228. Their basic characteristics are summarized in 
Table H 

Spectral types, apparent magnitudes, color indexes, 
parallaxes and proper motio ns we re all taken from the 
HIPPARCOS catalogue (ESA.il997l) . The bolometric cor- 



rections of all stars were derived from B — V hy using the 
calibration of Flower (1996 ) . 

The atmospheric parameters T^fr , logq, £t and [Fe/H ] 
and the stellar masses M* are from lSantos et all l)2003ll . 
The projected stellar rotational velocities, tisini, were 
obtained using the mea n ELODIE cross-corre lation func- 
tions (CCF) dip width ijBaranne et al."l996 V ex c ept fo r 
HD 106252, and the calibration bv iQueloz etal] l)l998|l . 
For HD 106252, vsmi has been derived fr om CCFs ob- 
tained with CORALIE and the calibration in lSantos et alJ 
12002} . 

Table |21 also gives for each star the equivalent width 
WA,Li of the A6707.8ALii line, as measured on the co- 
addition of all the available ELODIE spectra, corrected 
for the spectrograph straight light. It can be seen that all 
stars exhibit a lithium absorption feature which is rather 
strong except for two of them, HD 106252 and HD 190228. 
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HP 8574 



HP 23596 



HD 33636 



ELODIE 




3960 3965 3970 
X(A) 



3975 



3960 



3965 3970 
A{A) 



3975 



Fig. 5. The Can H line region for the six stars extracted 
from the co-addition of all spectra available for each star. 
The 3968.49 ACaii H line is indicated by a dotted Hnc. 
None of the stars exhibits any clear trace of chromosphcric 
emission feature in the line core (see text). 

Lithium abundances, scaled to log n(H) = 12, were deri ved 
following the method described in lNaef e t alJ ll2 001lJ) by 
using the curves of growth published bv ISoderb^fomet alJ 
l)l993(l and our tabulated Toff. The errors were estimated 
by assuming temperature and equivalent width variations 
of ± 1(7. The maximum abundance values are obtained by 
using Tcff -|- la and WA,Li + If and the minimum values 
correspond to T^s — Icr and W^xi — If - The resulting er- 
rors are not symmetrical because the error on Wa.li is 
supposed to be symmetrical and equal to ± 5 mA. 

The Can H absorption feature was extracted for the 
six stars as a chromospheric activity diagnostic. Fig.[3dis- 
plays the spectra of the Ca ii H line region obtained by the 
co-addition of all ELODIE spectra available for each star. 
None of the stars shows an appreciable emission. However, 
in this region, the ELODIE spectra do not usually have a 
high signal-to-noise ratio preventing us from the extrac- 
tion of meaningful values of the log{Rj^j^) indicator. The 
apparent absence of emission does therefore not preclude 
a mild activity. 

4. Orbital solutions 

The achieved precision (photon noise + instrumental er- 
ror) of the ELODIE radial-velocity measurements, al- 
though improved since the beginning of the ELODIE 
Planet Search Survey (see Sect. |21l, is of the order of 10 
ms~^ for these targets over the time span of the reported 
observations. The 254 individual high-precision ELODIE 



51600 52000 
JD-2 400 000 (d) 

Fig. 7. ELODIE and CORALIE radial-velocity mea- 
surements in kms^^ versus time for HD 33636. The 5 
CORALIE measurements are indicated by the arrow. Each 
measurement is plotted with an error bar equal to one es- 
timated standard deviation. Since the computed period 
is not yet covered, the orbital solution is not displayed 
for the whole fitted period. The orbital solution resulting 
from the best Keplerian fit, listed in TableOland discussed 
in the text, is sup erimposed (solid l ine). The orbital solu- 
tion published bvlVog t et al.l l|2002t) (dotted line) departs 
clearly from our solution at the epoch of the most recent 
measurements. 



radial velocities (-f 5 CORALIE measurements) obtained 
for the six planet-candidate host stars are available in elec- 
tronic form at the CDS in Strasbourg. All these measure- 
ments are displayed in Fig. Eland which also show the 
orbital solutions, resulting from the best Keplerian fit and 
described below for each star. 

We checked all solutions against false-alarms with 
Fourier transform analysis and found negligible (below the 
1% level) false-alarm probabilities. We also checked their 
fit errors with Monte-Carlo simul ations provided withi n 
the ORBIT programme described inl Forveille et aljl|l999tl . 
To this end, for each orbital solution, we ran 10 000 trials 
to derive Monte-Carlo estimates of the 1 and 3 a confi- 
dence intervals. None of the solution displays a noticeable 
discrepancy with these simulations except in three cases 
for a few parameters, as mentioned in the text below. 

4.1. HD8574 

We obtained a total of 41 ELODIE high-precision radial- 
velocity measurements for HD 8574 which was observed 
since Jan. 11, 1998 (HJD = 2 450 825). We list in Table |3 
the fitted orbital elements to these measurements together 
with the computed minimum mass 7712 sin i = 2.11 Mj^p 
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Fig. 6. ELODIE radial-velocity measurements in kms^^ versus time for all stars but HD 33636 shown in Fig. Each 
measurement is plotted with an error bar equal to one estimated standard deviation. The orbital solution resulting 
from the best Keplerian fit, listed in Table |2| and discussed in the text, is superimposed 



of the planetary companion, assuming a primary star 
mass of 1.17Mq. Our discovery of a companion of 
7712 sini = 1.95 Mj^p had been announced by the April 4*'* 
ESO Press Release^. The orbital solution presented here 
results from additional measurements obtained since then 
which further improved it. The errors listed in Table 01 are 
derived from Monte-Carlo simulations. They do not differ 



www. eso . org/outreach/press-rel/pr-2001/pr-07-01 .html 



appreciably from the fit errors assessing the quality of this 
solution despite the rather high value of Xred- 

The companion of HD8574 has the lowest period 
(227.55 d), well constrained (0.3%) and the best covered in 
terms of number of orbital cycles of the six planet candi- 
dates. Accordingly, it has also the smallest semi-major axis 
(a ~ 0.77 AU, i.e. 17mas). The slightly abnormal residuals 
from the fitted orbit yield a moderate "jitter" of 8.2ms~^, 
although the primary has no known activity and does not 
exhibit any clear one from the Can H line core (Fig.EJ. 
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Table 3. Fitted orbital elements to the radial-velocity measurements for the six stars, characteristics of their planet 
candidate companions and properties of the fit. When necessary, we replaced the fit error by the confidence interval 
value issued from Monte-Carlo simulations (see text). The m2 sini value is based on the listed in Table |21 ./Vmeas is 
the number of measurements that are effectively used in the orbital solution fit. (crv) is the mean measurement error 
and cr(0 — C) is the weighted rms of the residual. Comparison of the two values gives an indication of the amount of 
""jitter"" which can also be appreciated from the Xicd value (the reduced equal to x^/'^ where 1/ is the number of 
degrees of freedom) 
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5 Coralie measurements are included in the solution for HD 33636 (see text) 
^ at periastron 



at apoastron 



4.2. HD 23596 



of abnormal residuals and the agreement between fit errors 
and Montc-Carlo simulations. 



We obtained a total of 39 ELODIE high-precision radial- 
velocity measurements for HD 23596 which was observed 
since Jan. 10, 1998 (HJD = 2450824). Table 01 gives the 
obtained orbital elements together with the minimum 
mass of the companion (to2 sin i = 8.10 Mjup assuming 
Af* = 1.30Mq). The computed semi- major axis of the 
planetary orbit is a = 2.88AU, i.e. 55 mas. This solution 
seems to be well constrained as judged from the absence 



The planet of HD 23596 has the largest semi-major axis 
of the six candidate planets if we put aside the peculiar 
case of HD 33636. With a period of 4.28 yr, the candidate 
around HD 23596 belongs to the category of long period 
planets which is naturally growing with the time span of 
the on-going surveys. It is also a massive planet. We em- 
phasize the metal-rich nature of the parent star, a charac- 
teristics not shared by any of the six other stars. 
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4.3. HD 33636 

The fitted orbital elements for HD 33636 to its 42 ELODIE 
radial-velocity measurements obtained since Feb. 18, 1998 
(HJD = 2 450 863), plus 5 recent CORALIE measurements, 
are presented in Table El The velocity offset between the 
two instruments is an additional free parameter. Its fitted 
value is AiJ^E-c = 0.038 ±0.004 km s~^ The variability 
of the HD 33636 radial velocities has been evident for sev- 
eral years since the first measurements were done while 
the planet was at it s periastron. With a primary mass 
of 1.12 Mq given bv ISantos et al.l l)2003|) . we compute a 
minimum mass of m2 sini = 10.58 Mjup and a semi-major 
axis of a = 4.08AU, i.e. 142 mas, for the planet. It must 
be emphasized that the period is certainly not yet cov- 
ered (see Fig. CJ, a fact corroborated by the change of 
the fitted period with the addition of the measurements 
of our last observing campaign. At this stage, most of the 
constraints come from the impact of the eccentricity on 
the radial velocity curve shape. Therefore the period of 
7.74 yr is presently badly constrained but other param- 
eters should be less affected. Interestingly, Monte-Carlo 
simulations show that the fit errors (listed in Table 3) 
are somewhat overestimated. We adopt the asymmetric 
Monte-Carlo errors for the period and keep the fit errors 
for the other parameters. 

According to our preliminary orbital solution, 
HD 33636 b has the largest semi- major axis (4.08 AU) of 
the six planets. Despite its uncertainty, the period is long 
and therefore this new planet also belongs to the grow- 
ing category of long period planets. It is also be the most 
m assive of our six ca ndidates. 

IVogt et 311 112002^1 have announced the independent de- 
tection of a companion of m2 sini = 7.7Mjup. We have su- 
perimposed their orbital solution to ours in Fig. [7| The 
discrepancy between the two solutions is clearly due to 
their largely underestimated period. Our last campaign 
provides several measurements of good quality with which 
their orbital curve departs by 6 to 8 standard deviations 
of our most re cent measurements. Moreover, the period 
IVoet et al] |2002|) is shorter than the time span of our 



m ivoet et aij 11200 

measurements 



4.4. HD 50554 

A total of 41 ELODIE precise radial-velocity measure- 
ments were obtained for HD 50554 since Nov. 16, 1997 
(HJD = 2 450 769). We list in Table El the fitted orbital 
elements together with the computed minimum mass 
7712 sin i = 5.16 Mjup, assuming a primary star mass of 
1.11 Mq. The confidence interval resulting from Monte- 
Carlo simulations is a bit larger than the fit error for 
the period. This is probably due to the peculiar weight 
of the first season measurement in the period estimation 
(see Fig. jSJ . We therefore took the confidence interval for 
the period instead of its fit error. The semi-major axis 
of 2.41 AU, i.e. 77 mas, is derived for this new planetary 
companion. The mean radial-velocity measurement er- 



ror (eRv) and the weighted rms radial- velocity deviation 
(t(0 — C) differ notably, yielding an appreciable "jitter" 
of 6.3 ms^^. 

Our discovery of a companion of 7712 sini = 4.9 Mj^p 
had been announced by the ESO Press Release nb 07/01 
(April 4, 2001). The orbital solution presented here results 
from additional measurements obtained since then which 
fur ther improved it^ 

iFischer et al ] (120021^ have confirmed the detection of a 
companion with 7772 sin z = 3.7 Mjup. Except for their Ki 
and Lu values that differ significantly from ours, the two 
orbital solutions are compatible. The better coverage of 
the radial-velocity maximum by our measurements prob- 
ably explains the difference. The 77i2sin7 discrepancy re- 
sults from the rather large Ki difference between the two 

solutions^ 

The iFischer et al.l ((2002ll value of the chromospheric 
activity tracer log{Rjjj^) (—4.94), indicating a low activ- 
ity for HD 50554, is not useful for explaining the observed 
"jitter" , specially since the level of "jitter" is not well cor- 
related to the log (j?H-y) activity indicator f or low-activity 



relate d to the lom/tH-y I activity mdicator h 
stars l|Saar et alJll998[ ISantos et aLll2000b|) 



4.5. HD 106252 

We gathered a total of 40 ELODIE high-precision radial- 
velocity measurements for HD 106252 which was observed 
since March 1, 1997 (HJD = 2 450 509). In TablejSare pre- 
sented the fitted orbital elements to our measurements. 
With a primary mass of =1.02 Mq , we compute a min- 
imum mass of 77i2 sin 7 = 7.56 Mjup and a semi-major axis 
of 2.70 AU, i.e. 72mas, for HD 106252b. This solution is 
well constrained as confirmed by the absence of abnor- 
mal residuals. The Monte-Carlo confidence intervals are 
in good agreement with the fit errors. 

The discovery of a companion with 7712 sin 7 = 6.8 Mj^p 
had also been announced by the ESO Press Release nb 
07/01 (April 4, 2001). The orbital solution presented here 
results from additional measurements obtained since then 
which further improved it. With a period of 4.38 yr, this 
new planet also belongs to the growing category of long 
pe riod planets. It is als o a rather massive candidate. 

iFischer et ahi l|2002|) have confirmed the detection of 
a companion with 77x2 sin z = 6.96 Mjup, a value in rather 
good agreement with ours. Their period is at more than 
3(7 from our value. This difference probably results in their 
poorer temporal coverage. 

The Fisch er et ahl l)2002i) value of the chromospheric 
tracer log(i?^;^) (—4.97), indicating a low activity for 
HD 106252, agrees with the absence of abnormal residu- 
als. 



4.6. HD 190228 

We obtained a total of 51 ELODIE high-precision radial- 
velocity measurements for HD 190228 which was observed 
since Jul. 22, 1997 (HJD = 2 450 652). We list in Table |S1 
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Fig. 8. 1712 s'mi distribution for the 115 planet candidates 
discovered so far with m2 sinz < ISMjup. The grey square 
areas point to the contribution of the additional six new 
candidates presented in this paper 



the fitted orbital elements to these measurements together 
with the computed minimum mass m2 sin i = 3.58 Mjup 
assuming th e prirn ary star mass of 0.83 Mq given by 
ISantos et al.l l)2003j) . The semi-major axis a = 2.02AU, 
i.e. 33 mas, is given as well. The confidence interval re- 
sulting from Monte-Carlo simulations is somewhat larger 
(indeed asymmetrically) than the fit error for the eccen- 
tricity and Ki. This could be caused by the bad coverage 
of the periastron (see Fig. 0) . We therefore took the con- 
fidence interval for those two parameters instead of their 
fit error. 

The discovery of a co mpanion with m?. s in i = 3.4 Mjup 
had been announced in ISivan et all (|2000|). The orbital 
solution presented here results from additional measure- 
ments obtained since then which slightly improved it. 

Among the previous discoveries, more than a dozen 
planet candidates orbit sub-giant stars i.e. almost 20% 
of the total. Not yet statistically meaningful, this ratio 
however appears to be much larger than the proportion 
of sub-giant stars in the survey sample (19 stars among 
the ELODIE Planet Search Survey sample i.e. 5% of the 
total). The planet of HD 190228 is therefore not peculiar 
in this respect. On the contrary, HD 190228 is peculiar 
in terms of metallicity since it is the unique case in this 
paper sample with a metallicity value (—0.24 ± 0.06) quite 
below the solar abundance. It is also peculiar with respect 
to its spectral sub-type, since it is the coolest star of this 
sample. 



5. Discussion 

5.1. Existence of the companions 

The orbital periods derived for the six planetary compan- 
ions, spanning from 0.62 to 4.38 yr for five of them and 



even more for HD 33636 whose period is fairly imprecise, 
are much longer than the possible rotation periods of the 
primary stars, even though we have no individual estima- 
tion for most of them. This makes it very unlikely that any 
of the radial-velocity curves might be produced by activity 
spot phenomena as it h as been observed for the short pe- 
riod variable HD 166435 ()Queloz et al.l200lj) . Nevertheless 
we checked this for each star by computing the line bisec- 
tors of the observed CCFs and the linear correlation be- 
tween the radial-velocity measurements and the bisector 
span. With maximum absolute values of the correlation 
coefiicient below 0.23, we found no significant correlation, 
down to the la level, for any of the six stars. The interpre- 
tation of the Doppler signature by the presence of a com- 
panion thus casts no doubt for the six planet candidates, 
even for the mildly active stars, a conclusion in agreement 
with that of iNaef et al.1 l)2000aj) about long-period planet 
detectability around active stars. 



5.2. m2smi distribution and mass function 

With six new planet candidates, we can usefully up- 
date the distribution of m2sinz. Fig. [S] displays the his- 
togram of 7712 sin i of the 115 objects known so far, with 
7772 sin i < 18 Mjup, including our six new candidates and 
about thirty other new candidates announced very re- 
cently. These additional objects, half of which below 
7)7,2 sin i — 2 Mjup globally tend to smooth the cumulative 
frequency distribution and to fill s ome gaps present abo ut 
7)7,2 sm 7 — 8Mjup (see Fig. 1 of I.Torissen et al.l l(200lh . 
Because the techniqu es used to estimate t he deconvolved 
M2 distribution (e.g. Ijorissen et ahllioOlh is sensitive to 
the sampling quality, one can infer that such an approach 
will be even more robust with the new, larger input fre- 
quency distribution. More specifically, the latter is likely 
to result in a smoother M2 distribution, specially around 
and above lOMjup that could make it less obvious that 
10 Mjup would be a quasi upper M2 limit. Since this is only 
qualitative, it will be very interesting to know whether this 
study, applied to the new extra-solar planet sample, would 
rather lead to a vanishing distribution in the 15 Mjup re- 
gion. 



5.3. Eccentricity versus period diagram 

As seen above, the periods of all the companions are long, 
exceeding three years except for HD 8574. Their orbital 
eccentricities are close to 0.5 for four of the objects and 
about 0.3 for the two other ones. We have placed these 
six companions on a e vs P diagram together with all the 
other known companions with 7)7,2 sin 7 < 18Mjup (Fig.EJ. 
Their location is dictated by the large values of the pe- 
riods which is the natural effect of the bias towards long 
periods while the time span of the observing programmes 
increases. Apart from this, the diagram shows no particu- 
lar trend for the six companions which are spread in a re- 
gion of intermediate eccentricities, already well populated. 
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Fig. 9. The location of the 115 planet candidates discov- 
ered so far with TO2 sini < 18 Mjup in the diagram eccen- 
tricity versus period. The six new candidates are shown 
as open triangles 



The locations of our six new companions are in agreement 
with, and reinforce, the general behavior that tends to re- 
produce the spectroscopic binaries e vs P distribution for 
long periods l)Mavor et alJl200lHMazeh fc Zuckeil200l|) . a 
rather paradoxical constraint on planet formation mecha- 
nisms. 

5.4. Metallicity 

It has been shown that the metallicity distribution of the 
planet-bearing stars is statistically 0.25 dex higher tha. n 
for solar neighborhood stars (e.g. ISantos et alJ 
Since four of the six stars with new companions discussed 
here have solar metallicity, they will not help to statisti- 
cally improve this conclusion. Moreover, while HD 23596 
is over-metallic by 0.32 dex and therefore clearly falls in 
the upper part of the planet-bearing star distribution of 
ISantos et al1l)200l[l . HD 190228 is under-metallic by a sim- 
ilar amount (—0.24 dex). Both stars therefore almost com- 
pensates each other's effect on the distribution. The six 
stars will rather serve as additional material for further 
studies of metallicity distribution based on larger statisti- 
cal samples. 

5.5. Direct imaging feasibility 

With separations at the orbital elongation larger than 
60 mas, four out of the six new planet candidates could in 
principle be targets for diffraction-limited, i.e. adaptive 
optics, K band near infrared imaging on 8 m-class tele- 
scopes. The feasibility of such measurements is of course 
strongly dependent on the achievable contrast and there- 
fore very questionable for long periods, cool Jupiter-like 
planets. Such imaging would anyway place useful con- 
straints on the maximum mass and the orbital inclination 
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Fig. A.l. The ELODIE radial-velocity curve for 
HD 13507 versus time. The fit by the reflex motion due 
to the companion is superimposed to the measurements. 

of the companions. Altogether, the objective to image such 
companions clearly c alls for very high con trast dedicated 
instruments (see e.g. iMouillet et alJl200ll) . 
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Appendix A: HD 13507 

According to the HIPPARCOS catalogue, HD 13507 
is a my = 7.12 GO dwarf star with a parallax of 
38.12 ± 0.89 mas but its absolute magnitude Mv=5.10 
B -V (0.672 ±0.007) rather suggest a G4 spectral type. 

The atmospheric parameters and the mass of 
this star have been derived from a high signal- 
to-noise ELODIE spectrum by N.C. Santo s (priv . 
comm.) in the way described in ISantos et al.l l)2000aj) : 
[Fe/H] = 0.01 ± 0.08, T^s = 5775 ± 80 K, log g = 4.71 ±0.17 
(cgs), ^t = l-33±0.11kms-i and = 1.09 Mq. 
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A total of 19 ELODIE precise radial-velocity mea- 
surements were gathered for HD 13507 since Jan. 8, 1998 
(HJD = 2 450822). A preliminary but apparently good so- 
lution, based on a linear drift of 341 ± 8 ms~^ yr~^ and a 
Keplerian variation compatible with a planetary compan- 
ion, that was announced as such, could be derived until 
mid 2002. 

The later measurements obtained between July and de- 
cember 2002 (see Fig. lA.ljl invalidated this interpretation 
and revealed instead a classical spectroscopic binary veloc- 
ity curve, caused by a low mass star companion. This pos- 
sibly explains the noticed discrepancy between the spec- 
tral type from HIPPARCOS and its photometry. The com- 
puted orbital parameters obtained from a still preliminary 
fit are: period PciSOOOdays, e = 0.14, ivTi = 0.694 ms-\ 
minimum mass of the companion m2 sini = 52 Mjup, as- 
suming a primary star mass of 1.09Mq, and semi-major 
axis a = 4.3AU, i.e. 164 mas. This solution is superim- 
posed to the measurements, displayed in Fig. lA.ll 

We have searched for a companion to HD 13507 by 
adaptive optics imaging. An image was obtained in 
Bracket7 with PUE'O at the CFHT on August, 7*'' 2001 
(HJD = 2 452 128.61) but no companion shows up at a 3(T 
level. A unique negative measurement leaves place to a 
possible alignment of the components but monitoring the 
binary with high dynamics imaging over a few years should 
resolve it and permit to derive the companion mass, a use- 
ful information for studies of sub-stellar objects statistics. 
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